Gamma-ray emission from strongly magnetized pulsars 
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ABSTRACT 

In a strong magnetic field, B > 4 x 10 12 G, 7-rays emitted nearly along curved 
field lines adiabatically convert into bound electron-positron pairs (positronium 
atoms) rather that decaying into free pairs. This process may modify the polar 
gaps of strong magnetized pulsars. Unlike free pairs, such bound pairs do not 
screen the electric field component along the magnetic field in the polar gaps. As 
a result the total power carried away by both relativistic particles and radiation, 
from the polar gap into the pulsar magnetosphere, may increase significantly 
(up to a few tens times) in comparison with the conventional polar gap models 
where creation of bound pairs is ignored, and it may be a substantial fraction 
of the spin-down power. We demonstrate that the total power of the modified 
polar gaps may be enough to explain the observed non-thermal luminosities of 
all known strongly magnetized, 7-ray pulsars. 

Subject headings: gamma rays: theory — pulsars: general — radiation mecha- 
nisms: non-thermal 



1. Introduction 



After the Compton Gamma Ray Observatory (CGRO) ended its activity in 2000, seven 
rotation-powered 7-ray pulsars were known (e.g., Thompson et al. 1997). More recently, 
the AGILE telescope reported the detection of two other radio pulsars in 7-rays. The Fermi 
Gamma-ray Space Telescope was successfully launched in 2008. After two years of observa- 
tions by Fermi, the number of know 7-ray pulsars increased dramatically, to more than sixty, 
and continue increase (Ray & Saz-Parkinson 2010). For known 7-ray pulsars the maximum 
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of the radiated power is concentrated in the 7-ray range, except for the Crab pulsar and PSR 
B 1509-58, the power of which peaks in hard X-rays. The pulsar luminosities L 7 in 7-rays are 
a substantial fraction of the spin-down power, rj^ = L 1 /E TOt ~ 10~ 2 — 1 (Abdo et al. 2010), 
where E rot = 4tt 2 IPP~ 3 is the spin-down power, / is the moment of inertia of the neutron 
star (generally taken to be 10 45 gcm 2 ), P its spin period, and P its spin-down rate. 

The starting point common to all viable models of pulsars is that strong electric fields 
are generated in the magnetospheres of rotating neutron stars (e.g., Michel 1991). The 
electric field component, En = (E-B)/|B|, along the magnetic field, B, is non-zero and 
may accelerate particles to ultrarelativistic energies. The accelerated particles emit 7-rays 
due to curvature emission and other processes. Some of the 7-rays are converted into e + e~ 
pairs in a strong magnetic field. The pairs screen the field E\\ in the pulsar magnetosphere 
everywhere except for some compact regions called "gaps". These gaps are, in fact, "engines" 
responsible for the non-thermal radiation of pulsars. 

Three kinds of gaps have been proposed in trying to explain the non-thermal radiation 
of pulsars. The main difference between these is in the site of gaps. A gap that forms near 
the magnetic pole of a pulsar is called a polar gap (Ruderman & Sutherland 1975). Besides, 
an outer gap may form between the surface of null Goldreich- Julian density and the light 
cylinder (Cheng, Ho, & Ruderman 1986a,b). The third possibility is a slot gap located in the 
space-charge-limited flow along the last closed field line (Arons 1983; Muslimov & Harding 
2004). 

In conventional models of polar gaps, it is assumed that created e + e~ pairs are free 
(Ruderman & Sutherland 1975; Daugherty & Harding 1982, 1996; Medin & Lai 2010 and 
references therein). However, this assumption is not valid if the magnetic field is higher than 
~ 0.1 B Q ~ 4 x 10 12 G, where B = m 2 c 3 /eh = 4.4 x 10 13 G. The reason is that in such 
a strong field, 7-rays emitted nearly along the curved field lines, adiabatically convert into 
bound electron-positron pairs (positronium) rather than decaying into free pairs (Shabad & 
Usov 1985, 1986; Herold et al. 1985; Usov & Melrose 1995, 1996; Baring & Harding 2001; 
Harding & Lai 2006; Thompson 2008). The fact that e + e~ pairs created in a magnetic field 
B > 0.1B are bound, may be very important for the observational appearances of strongly 
magnetized pulsars (Usov & Melrose 1995). In particular, unlike free pairs, such bound pairs 
do not screen the electric field E\\ near the pulsar, which requires a net charge density to 
build up. As a result the pulsar luminosity is higher than it would have been if the created 
pairs are free. For pulsars with strong surface magnetic fields, B s > 0.1B CT , a polar gap is 
modified by taking into account creation of bound pairs as was suggested by Usov & Melrose 
(1995, 1996). In this model, the total power carried away by both relativistic particles and 
radiation, from the polar gap into the magnetosphere, may increase significantly (up to a 
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few tens times) in comparison with the conventional polar gap models, and it may be a 
substantial fraction of the spin-down power. In this Letter, we demonstrate that the total 
power of the modified polar gaps may be enough to explain the observed luminosities of all 
known strongly magnetized, 7-ray pulsars. Other observational consequences of the modified 
polar gap model for these pulsars are also discussed. 



2. Radiative efficiency 

In the modified polar gap model (Usov & Melrose 1985, 1986), the fraction of the spin- 
down power carried away by both relativistic particles and radiation from the polar gap into 
the magnetosphere is estimated as 

*-£-!«)-[-«)■ 

where Pi ~ 0.5(P P /0.1 Po) 2 ^ 3 s j an d P p is the magnetic field strength at the poles. At 
P = 2~ 2 / 3 P 1 ~ 0.63 Pi, the value of 77^ is maximum, rj^\ max = 3/8 = 0.375. In this case the 
polar gap luminosity L b p is comparable with the spin-down power. 

The modified polar gap model is valid if both P p > 0.1B and P 2 < P < P\ (Usov & 
Melrose 1986), where 

P 2 ~ 0.07(T S /10 6 K) 4 / n (P p /0.1P cr ) 2/11 s. (2) 

For 7-ray pulsars satisfied these conditions, measured and derived properties are given in 
Table I. Additionally, the parameters of two 7-ray pulsars (J0633+1746 and J1057-5226) are 
included in Table 1. In these two cases, the observed efficiency of transformation of the spin- 
down power into the non-thermal high-frequency emission is very high, ?7° bs = L 1+ x / E rot ~ 
1, while the dipolar B-field estimate is only slightly below the value required for bound-pair 
formation. The surface magnetic field may be higher than the dipolar estimate and satisfies 
the requirement of bound-pair formation, B p > O.IPoj provided one invokes higher-order 
multipolar components, or an off-centered dipole. To estimate the possible value of 77^ for 
these two pulsars we take B p = O.IBq = 4.4 x 10 12 G (see the second lines in Table 1). 

From Table 1 we can see that the predicted value of 77^ is about or more than the 
minimum observed value of the radiative efficiency, 77° bs , for all strongly magnetized 7-ray 
pulsars except of PSR J2021+4026. It is worth to note that the 7-ray luminosities of pulsars 
given in Table 1 were only roughly estimated by Abdo et. al. (2010) using the following 
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equation 



L, = ^d 2 f n G, 



(3) 



where d is the distance of the pulsar, G is the average energy flux in 7-rays, and /q is the 
flux correction factor that takes into account the angular anisotropy of the 7-ray emission of 
pulsars. The value of /q is model-dependent and varies from ~ 1 in the outer gap and slot 
gap models to ~ l/4ir ~ 0.08 in the polar gap model (e.g., Thompson et al. 1994). Abdo 
et al. (2010) have used fn = l throughout the paper, and therefore, the 7-ray luminosities 
may be overestimated at least a few times. Summarizing, the power carried away by both 
relativistic particles and radiation from the modified polar gaps into the magnetosphere may 
be sufficient to explain the non-thermal luminosities of all known strongly magnetized 7-ray 
pulsars, including PSR J2021+4026. 

For comparison, in conventional polar-gap models, where the created pairs are assumed 
to be free, the corresponding fraction of the spin-down power going into both high-energy 
particles and radiation may be estimated as (e.g., Usov & Melrose 1995) 



From Table 1 we can see that rjf is more than an order of magnitude smaller than r}° hs 
except for a few 7-ray pulsars. Even taking into account the uncertainty in i]° hs because of 
the uncertainty of fa, we can see that for the main part of the strongly magnetized 7-ray 
pulsars the inferred high efficiency of conversion of rotational energy into 7-ray radiation 
cannot be explained within the conventional polar-gap models. 

For six pulsars the observed luminosities in 7-rays are not presented in Table 1 because 
their distances are unknown. Equalizing ?7° bs to 77^ for these pulsars, upper limits on the 
distances may be estimates in the frame of the modified polar gap model. 



The 7-ray spectra of pulsars were fitted by an exponentially cutoff power-law model, 
dN/dE (x E~ a exp(—E / E cuoS ), with the photon index a in the range 0.7-2.4 and the cutoff 
energy E cuoS at 0.1-7 GeV (Abdo et al. 2010, Pilia et al. 2010). These spectra may be 
explained by curvature radiation of high-energy electrons that escape from the polar caps 
with the Lorentz factor T ~ 10 7 (Harding 1981, Daugherty & Harding 1982; Chiang & 
Romani 1992). The value of T ~ 10 7 is more or less appropriate for primary electrons at 
the top edge of the polar gap in the conventional polar gap models, in which the density 





3. Spectrum of 7-rays 
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of primary particles is slightly smaller than the Goldreich- Julian density n G j (Goldreich & 
Julian 1969; Ruderman & Sutherland 1975; Arons 1981). To explain the luminosities of 7-ray 
pulsars, it was suggested by Harding (1981) that the density of the outflowing electrons with 
Lorentz factor T ~ 10 7 is ~ (10 2 — 10 3 )ncj. Although this suggestion is inconsistent with 
the conventional polar gap models, it may be consistent with the modified polar gap model 
(Usov & Melrose 1995, 1996). Indeed, inside the modified polar gaps the primary electrons 
are accelerated by the electric field E\\ ~ E^° n to the Lorentz factor r prim ~ (2 — 3) x 10 s and 
generate very hard 7-rays via curvature radiation, where E™ n ~ (1 — 2) x 10 10 V/cm is the 
electric field at which ield ionization of bound pairs becomes important. These 7-rays are 
absorbed in the magnetic field and create secondary electron-positron pairs with the mean 
Lorentz factor and the density 



respectively (Usov & Melrose 1996). The parameters of the secondary particles are well 
consistent with those suggested by Harding (1981) to explain both the luminosities and 
spectra of 7-ray pulsars. 

In the frame of the conventional polar gap models, it was argued that photon splitting, 
7 + B — > 7' + 7" + B, and magnetic absorption of photons, 7 + £? — > e + + e~ + B, have to 
reduce the cutoff energy E cutoS for strongly magnetized 7-ray pulsars, B p > 10 13 G (Harding 
et al. 1997; Baring & Harding 2001; Baring 2004). This is consistent with the observational 
data on PSR B1509-58 that has one of the highest magnetic fields, 5 P ~ 3.1 x 10 13 G, 
and the softest spectrum with a cutoff at E cuto s = 81 ± 20 MeV. The second softest 7- 
ray spectrum with E cuto s = 0.7 ± 0.5 GeV is that of PSR B0656+14, recently observed by 
Fermi. This pulsar has a rather strong magnetic field, B p ~ 0.93 x 10 13 G, and its soft 7-ray 
spectrum may be explained in the same way. At present, however, there are a handful of 
7-ray pulsars with surface magnetic fields higher than that at the surface of PSR B0656+14, 
while their 7-ray emission is present to at least a few GeV. One such pulsar is LAT PSR 
J0007+7303, in the supernova remnant CTA 1. This pulsar has one of the hardest 7-ray 
spectra with E CVLtoS = 4.6 ± 0.4 GeV, while the magnetic field at its poles is B p ~ 2.2 x 10 13 
G, which is only slightly smaller than the _B p -value for PSR B1509-58. Among other pulsars 
with hard (E cnto s > a few GeV) 7-ray spectra and strong (B p > O.IBq) magnetic fields 
are PSR J0631+1036, J0633+1746, J1709-4429, and J2021+4026. Since one expects strong 
absorption of hard 7-rays in the vicinity of these pulsars (Harding et al. 1997, Baring 2004), 
it is concluded that (1) a polar gap model is inconsistent with the available data on the 
strongly magnetized 7-ray pulsars with hard 7-ray spectra, and (2) the 7-ray emission of 
these pulsars arises largely in the outer magnetosphere (e.g., Abdo et al. 2010; Michelson, 
Atwood, & Ritz 2010). 
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However, these conclusions may be premature at least for some strongly magnetized 
7-pulsars. Indeed, in the modified polar gap model 7-rays generated near the surface of a 
strongly magnetized pulsar adiabatically convert into bound e + e~ pairs. If the photoioniza- 
tion of the bound pairs is small, these pairs mostly flow away from the pulsar vicinity and 
may annihilate at large distances generating hard 7-rays that escape from the pulsar mag- 
netosphere. In other words, bound pairs may be an intermediate agent that transfers 7-rays 
from the pulsar vicinity to the outer magnetosphere where their absorption is negligible. 

The surface temperature of PSR B1509-58 may be as high as ~ 2 x 10 6 K. Besides, 
this pulsar is a powerful (Lx — 10 35 erg/s) source of soft X-rays with a non-thermal (power- 
law) spectrum. The soft X-ray and 7-ray pulses are phase-aligned, indicating that they 
are generated in the same region in the pulsar magnetosphere (near the magnetic poles in 
the polar gap model). Most probably, the bound pairs are mostly ionized by both thermal 
and non-thermal X-rays before leaving the vicinity of PSR B1509-58. In this case, the 
formation of bound pairs does not affect the conclusion of Harding et al. (1997) that in 
the polar gap model the 7-ray spectrum of PSR B1509-58 is expected to be very soft. As 
to PSR J0007+7303, there is an upper limit of T s < 6.6 x 10 5 K (Halpern et al. 2004). 
Recently, Caraveo et al. (2010) decreased the upper limit on the surface temperature of PSR 
J0007+7303: T s < 5.3 x 10 5 K, that is one of the most constraining data points on cooling 
models of neutron stars. This makes PSR J0007+7303 by far the coldest neutron star for 
its age interval, suggesting the necessity of enhanced neutrino emission for this rather young 
7-ray pulsar (Page, Geppert, & Weber 2006; Caraveo et al. 2010). The surface temperature 
of PSR J0007+7303 in the cooling model of neutron stars with enhanced neutrino emission 
may be as low as (1 — 2) x 10 5 K (Page et al. 2006). For this temperature, the bound pairs 
may leave the pulsar vicinity and annihilate far from the pulsar into hard7-rays, £ 7 > a few 
GeV. 



4. Discussion 

In this paper, we have demonstrated that the observed non-thermal luminosities of 
strongly magnetized pulsars may be explained in the polar gap model modified by taking 
into account the process of adiabatical conversion of 7-rays into bound electron-positron 
pairs in a strong magnetic field, B > 4 x 10 12 G. The high-energy spectra of 7-ray pulsars 
may be also explained by the modified polar gap model. However, all available polar gap 
models have difficulty in reconciling the observed 7-ray pulses with the known radio pulsar 
geometry. At present, it is not clear possible or not to overcome this difficulty in the modified 
polar gap model at least for some strongly magnetized pulsars. 
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Table 1: Properties of 7-ray pulsars 
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